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Abstract 

We investigate the gauged Nambu-Jona-Lasinio model in curved spacetime at the large N c limit 
and in slow-roll approximation. The model can be described by the renormalization group cor¬ 
rected gauge-Higgs-Yukawa theory with the corresponding compositeness conditions. Evaluating the 
renormalization group (RG) improved effective action, we show that such model can produce CMB 
fluctuations and find inflationary parameters: spectral index, tensor-to-scalar-ratio and running of 
the spectral index. We demonstrate that the model can naturally satisfy the Planck 2015 data and 
maybe considered as an alternative candidate for Higgs inflation. 
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1 Introduction 

The standard model of the particle physics is based on the gauge symmetry and the symmetry breaking. 
For example, QCD is described as SU(3 ) gauge theory. The gauge interaction causes chiral symmetry 
breaking. Because of the strong coupling constant, low energy phenomena are often approximated by 
using effective models. For instance, Nambu-Jona-Lasinio (NJL) model [T] is a well-known model for 
light mesons. Gauged NJL model maybe applied to the description of QCD, as SU{N C ) gauge theory, 
in which the chiral symmetry is broken and meson like composite fields are generated at high energy 
scale (see, for example [2|)- The renormalizability of the model is precisely studied by using the ladder 
Scliwinger-Dyson equation in flat spacetime [3]. In the leading order of the 1/N C expansion the model 
can be described by a gauge-Higgs-Yukawa type Lagrangian with compositeness conditions [4] . 

Gauged NJL model is expected to be very important at the early universcQ. The vacuum energy for a 
fermion field is negative, unlike to the scalar field. After the dynamical symmetry breaking a composite 
scalar can develop a positive vacuum energy. A primordial potential energy for the composite scalar may 
induce inflation, i.e. exponential expansion of our universe. An evidence of inflation is observed in the 
cosmic microwave background radiation (CMB) [6, 2- In this paper we propose gauged NJL inflation 
model as an alternative to Higgs inflation. We evaluate the CMB fluctuations in the gauged NJL model 
and demonstrate that inflationary parameters are consistent with Planck bounds. 

The study of NJL model in curved space-time has been started some time ago (for review, see Ref. 0). 
The curvature coupling for the composite scalar field has been obtained in Ref. mm- In the gauged NJL 
model it has been shown that the chiral symmetry can be broken for a strong curvature m- In Ref. [12] 
the authors have discussed the contributions from the dynamical symmetry breaking to the expansion 

1 Constraints from the high energy collider physics are reviewed in Ref. [5]. 
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of the universe. However, as far as we know there is only preliminary work to argue a contribution 
from composite fields to the CMB fluctuations. For instance, in Ref. jT3] the authors have assumed 
non-vanishing fermion condensation and evaluated the CMB fluctuations. 

This article is organized as follows: In Sec. [2] the gauged NJL model is introduced. The model is 
rewritten as renormalization-group corrected gauge-Higgs-Yukawa Lagrangian with compositeness con¬ 
ditions. In this way, this looks qualitatively similar to the study of RG improved effective potential in 
curved spacetime 2] and corresponding RG improved inflation [T5]- In Sec. [3] we calculate the infla¬ 
tionary parameters in the gauged NJL model under the slow -roll approximation scenario. In Sec. 0] the 
inflationary parameters are numerically calculated and compared with the Planck data. Then we find a 
possible parameters range consistent with Planck data for the gauged NJL model. In Sec. [5] we briefly 
discuss an exit from inflation. Some concluding remarks are given in Sec. [6] 


2 Gauged NJL model 


We consider Nf fermion flavors with SU ( N c )<giQ gauge symmetry. It is assumed that the gauge interaction 
for Q is strong enough and well described by four-fermion interaction. We employ the gauged NJL model 
as an effective model of the SU(N C ) <8> Q gauge theory. The Lagrangian is given by 


-gNJL 


= C 


gauge 


■ -i pijfii!) + 


167 T 2 g 4 

8N f N c A 2 L 


+ (4>i'Y5T a 'tp)‘ 


( 2 . 1 ) 


where r“ are the generators for the flavor symmetry, C gaU ge shows the Lagrangian of the SU(N c ) pure 
gauge theory and D is the covariant derivative for the SU(N C ) gauge invariant and generally covariant 
fermion kinetic term. The gauge field for the gauge symmetry Q is integrated out. We introduce the 
compositeness scale A and use it to define the dimensionless four-fermion coupling g 4 . It is noted that 
the Lagrangian has the U{Nf)i, ® U{Nf)n flavor symmetry. Thus the model is regarded as a scale 
up of QCD with QED interaction. The fermion fields, ijj, are confined in baryon like states above the 
compositeness scale. We consider the model whose compositeness scale is higher than the one for inflation. 
The composite fields dynamically acquire compositeness scale mass. Only the possible environment to 
test the model is found at early universe. 

According to the auxiliary field method the Lagrangian (12.11) is rewritten as 


r — r 

t - / aux — a ~ , q 


■ ip [i.Ip 


cr — Z 7 5 t 7 r 
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Replacing the auxiliary fields cr and 7r a by the solutions of the equations of motion, 

167r2 94 7( a 16tt 2 54 

cr = -TTTTTT7#, 7T = - . 9 If), 


AN f N c A 2 


ANfN c A? 


( 2 . 2 ) 


(2.3) 


the original Lagrangian m is reproduced. Starting from the Lagrangian m at the compositeness 
scale A, one obtains the effective Lagrangian with a gauge-Higgs-Yukawa form below the compositeness 
scale through the renormalization group evolution [3], 


C g HY = + ^ad^a + \d^ a 8^ a - V + ttV) - ^(a 2 + ttV) 2 

—-^R(ct 2 + 7T a 7T a ) + ppilfilp — + i"/5T a TT a )lp- 


(2.4) 


Because of the U(N/)l <8> U(Nf)ii flavor symmetry one can assume that only the composite scalar field, 
cr, contributes the vacuum energy and induce exponential expansion of the universe. Below we drop the 
pseudo-scalar field, n a . 

The RG equations for the gauge coupling g(t), the Yukawa coupling y(t), the quartic scalar coupling 
A (t) and the curvature coupling £(t) are found in the leading order with respect to a modified 1/N C 
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expansion [TH[T 6 ] . 


16772 dt = 

1671-2 ^2/W = J/(0[«2/ 2 (i) - cg 2 (t)], 

16n2 ^i X ^ = 4 a 2 / 2 Wt A W ~ y 2 ( 01 > 


1671-2 = 2 ay 2 (i) 
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(2.5) 


where RG parameter t is defined by t = ln(yLi//Lio) with the renormalization scale /i and a reference scale 
/io, a, b, c and u being the positive constants, 

1 N 2 - 1 

a = 2 N f N c , b = -(11N C - 2N f ), c = 3^—. (2.6) 

In Ref. (IT it has been shown that the Lagrangian (12.41) is identified with ( 12 . 21 ) under the compositeness 
conditions in a fiat space-time, 


1 ^ n A(f A ) ^ m 2 (t A ) 2a 2 / 1 1 \ 

3 2 (Ay) ’ 3 4 (^a) ’ V 2 {tA) 1671-2 \34 tt(t A )) ’ 


(2.7) 


where fl(fA) is a function of t and it has a constant value lo at the limit, b —> +0 [TB]. It is shown that 
the composite scalar is conformally coupled to curvature [9j- Thus the compositeness condition for £ is 
given as, 

£(*a) = l- (2.8) 

Solving the RG equations (E2D with the compositeness conditions m and m, we obtain the 
running couplings, 


y 2 (t) - 
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y 4 (t) 
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where a(t ) is the running fine structure constant for the SU(N C ) gauge interaction, a{t) = g 2 (t)/( 47 t). 
At the fixed gauge coupling limit, b —> +0, Eqs. (12.91) and ( 12 . 101 ) reduce to 
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where w is defined by w = 1 — a/(2a c ) with a c = 2 t:/c. We also define y* and A» at the large A limit. 
At this limit we can explicitly solve the RG equation (12.51) for m 2 {t) and find 


mn 2 {t) = 


2 a 
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where gm(g) denotes the renormalized four-fermion coupling and g\ R is a constant parameter, see 
Ref. [TB]. Below we consider only the fixed gauge coupling limit, b = 0, for simplicity. 
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let us consider CMB fluctuations in the gauged NJL model. For this purpose, we evaluate the RG 
improved effective potential for the composite field ct in curved spacetime. We follow the procedure 
developed in Ref. HU. The one-loop effective potential for the theory (12.41) up to terms linear in R is 
given by 


V llo °P(a) 


—m 2 a 2 + -Act 4 + -£I?ct 2 
2 4 2 


ay 4 a 4 y 2 a 2 
2 • 16^ 2 ^ n ^ 


aRy 2 a 2 
12 • 16vr 2 
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(2.15) 


To obtain the RG invariant effective potential of the gauged NJL model one introduces the compositeness 
condition into the effective potential (12.151) through the RG equations. The solution of the RG equation 
for the effective potential satisfies 


V(g, y, A, to 2 , £, ct, y) = V ( g{t ), y(t), A(t), m, 2 (t), £(f), /xe*), 


(2.16) 


where the barred quantities are the renormalized ones at the scale y.e 1 . The scale is fixed to drop the 
logarithmic terms in the RG invariant effective potential [T5] , 


- ( f) 


l/(2-w) 


(2.17) 


It should be noted that the gravitational couplings are not represented in Eq. (12.161) . These terms do not 
contribute at the linear curvature approximation. 

Using the one-loop effective potential (12.151) as boundary condition, we calculate the RG invariant 
effective potential (12.161) from Eqs. (12.121) . (I2.13D . (12.141) and (12.171) . At the limit A —» oo the RG invariant 
effective potential is given by 
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(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


Thus, the gravitational effective action for the composite scalar field with weak curvature R <C ct 2 in the 
gauged NJL model is obtained, 


S = 


d A X\J~^g 


-V(a,n) 


( 2 . 22 ) 


where the reduced Planck mass is set as M p = (87rG) 4 / 2 = 1. 


3 Inflationary parameters in gauged NJL model 

The gauged NJL model may give us an alternative scenario to induce the exponential expansion of our 
universe. Here we assume that the typical scale of the inflation is much smaller than the compositeness 
scale, A, and consider the effective potential (12.181) for the composite scalar field, ct, as an inflaton 
potential. Then we calculate the inflationary parameters generated by the composite scalar field. 

For practical calculations it is more convenient to change the frame into the Einstein frame where the 
interaction term between the composite scalar and the curvature disappears. This change of the frame 
is achieved by the Weyl transformation, 


= n(x)~ 2 g^, 


(3.1) 


4 


































where g^ is the metric tensor in the transformed frame and the Weyl factor, f2(a0, is an analytic 
function with respect to the space-time coordinates. Applying the Weyl transformation (13.11) and setting 
the conformal factor to be 


n2 = i + ^)( CT 2 ) l/(2-«) } 


(3.2) 


one gets the effective action (12.221) in the Einstein frame 
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where g and R are the determinant of the metric tensor and the Ricci scalar in the transformed frame, 
respectively. In order to obtain the canonical kinetic term we redefine the auxiliary field cr to satisfy the 
relation, 


dip = 


2 \fl 2 2 — w 6 
Then the effective action m reduces to 


1 2 D{p) ^2\w/(i-2w) 
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1/2 


da 


(3.4) 


(3.5) 


Se = j d^xsf-'g - V E (a, g)^j , 

where the subscript E denotes the Einstein frame and the effective potential V E (a, g) is defined by 

V E (cr, g) = IT 4 [^a 2 + ^M( ct 2 ) 2/(2- w) 1 _ (3 . 6) 


The effective potential V E (cr, g) in the Einstein frame depends on the parameters in the original gauged 
NJL model through Eos. (12.191) - (12.211) . The quantum effects are encoded via the /z-dependence of 
these parameters with the compositeness condition. It should be noted that the obtained inflationary 
parameters are just the same in both the original and Einstein frames Ha- 

We will not be interesting in the exact solutions for the action (13.51) as it has been discussed in detail 
for various choices of the potential. We will concentrate here on quasi-de Sitter phase which eventually 
describes the inflationary dynamics. In the study of quasi-de Sitter phase we apply the standard slow- 
roll approximation which is well-known (see Refs. (19) . for instance) to be consistent approximation to 
describe the quasi-de Sitter phase. 

We consider a locally flat Friedman-Robertson-Walker (FRW) universe defined by the metric, 

ds 2 = dt 2 — a 2 [t)(dx 2 + dy 2 + dz 2 ), (3.7) 

where the evolution of the spacetime is described by the temporal development of the scale factor, a(t). 
In the FRW universe, the equation of motion for ip and the Freedman equation are given by 

dV E 

ip + 3Hip =——, (3.8) 


H 2 


1 fl 


3 \ 2^ 


V E 
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with the Hubble parameter, H = a/a, respectively. In the slow roll scenario we assume that the spacetime 
is slowly varying and impose the conditions 


ip 2 < V, ip <€. 


dV E 

dp 


(3.10) 
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Then we obtain the quasi-de Sitter phase, 


a(t + At) ~ xa(t)e'^ At + (1 — x)a(t)e V^T At ; 


(3.11) 


where x is a constant parameter. 

In a standard scenario of the the inflation the CMB fluctuation arises from the primordial quantum 
fluctuation for the inflaton field. Here we employ the slow-roll approximation for the inflationary expan¬ 
sion of the universe. Under the approximation the inflationary parameters can be fully represented by 
means of the inflaton potential, 
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The e-folding number N is given by 


da 
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N = 


Ve 
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V E 
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(3.13) 
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(3.15) 


where the integration is performed on the interval from the horizon crossing a = cx/v to the end of the 
inflation a = a en d . Since the inflationary parameters have to be small, the field configuration, a en d, can 
be estimated where one of the inflationary parameters exceeds order one. We noted that the change of 
the value for a en d gives only a small correction to the e-folding number N. The dominant contribution 
comes from the vale at the horizon crossing, cr/v, where the derivative of the effective potential is small 
enough. In the other words the field configuration, a n is obtained as a function of the e-folding number 
N from Eci. (I 3.15D . Then one can evaluate the density fluctuation, 6 , by [19] 

(3.16) 

<7 —a AT 


s ~ Ve 2 { dV E da \ ~ l 
\J\ 2 e 2 V da dp) 


The spectral index, n s , the tensor-to-scalar ratio, r, and the running of the spectral index, a s , are 
represented by using the slow-roll parameters ED], 


n s -l = 2 n | a=(rjv - 6 e| a=CTN , 
r = 16 el_, 


a s = 


dn s 
din k 


= -24 e z 


+ 16 eq\ — 2 £| 


a=a N 5 


(3.17) 

(3.18) 

(3.19) 


where we omit higher-order terms. Thus the observed CMB constraints for <5, n s , r and a s restrict the 
possible range of the parameters in the gauged NJL model for quasi-de Sitter universe under discussion.. 


4 Numerical results 

The effective action (12.221) depends on the fine-structure constant a, the renormalized four-fermion cou¬ 
pling < 74 /j, the number of fermion species, Nf, N c and the renormalization scale p,. In this section we 
numerically calculate the inflationary parameters 5, n s , r and a s with varying the model parameters, Nf, 
N c , g 4 E, a, p and the e-folding number, N. 

4.1 Small number of fermion species (Nf — 1, N c — 10) 

In this subsection we consider a model with a small number of fermion species, Nf = 1, N c = 10. First 
we set model parameters as l/g 4 R — I/sIr = 1> M = 1, and study the gauge coupling dependence. The 
behavior of the effective potential is drawn in Fig.|T]for a = 0.5 x 10 -12 ,1 x 10 -12 and 2 x 10 -12 . We plot 
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the points, a = a en d , a = &n =50 and a = o'n- 60 on each lines. A gentler slope of the potential is observed 
for a weaker gauge coupling. It should be noticed that we have to fine-tune the field configuration at the 
horizon crossing to obtain the enough e-foldings in the case of the NJL model (See Appendix A). The 
density fluctuation decreases as the slope becomes gentle. As is shown in Fig.( 2 j the density fluctuation is 
monotonically increasing as a function of a. The observed density fluctuation S ~ 4.93 x 10 -5 is obtained 
for a ~ O( 10 -12 ). 


3. x 10‘ 10 
2.5 x 1CT 10 
2. x 1CT 10 
w 1.5X1CT 10 
^ 1.x10“ 10 

5. x 10 -11 
0 

0 5 10 15 

<J 

Figure 1: Behavior of the potential Ve(ct,/ x) as a 
function of a for Nf = 1 , N c = 10, 1 / g^R— 1 / g\ R = 
1, ii = 1 and a = 0.5 x 10“ 12 ,1 x 10~ 12 ,2 x 10" 12 . 



a 

Figure 2: Behavior of the density fluctuation <5 as a 
function of a for Nf = 1, N c = 10, 1 /<? 4 .r— 1 / 54 # = 
1 and n = 1 . 



In Table |T] the values for the inflationary parameters n s , r and a s for a = 10 -12 ,10 -11 and 10~ 10 
are shown. The spectrum index n s , tensor-scalar ratio r, and the running of the spectral index, a s , are 
consistent with the observed Planck constraints, n s = 0.9655 ± 0.0062, a s = —0.0084 ± 0.0082 (68% CL, 
PlanckTT+lowP) and r 0 .oo 2 < 0.10 (95% CL, PlanckTT+lowP)[?/ It is also found that the a-dependence 
for the inflationary parameters n s , r and a s is negligible. 


a 

s 

n s 

r 

a s 

10~ 12 

3.58 x 10" 5 

0.961 

0.0084 

-0.00077 

10- 11 

1.13 x 10" 4 

0.961 

0.0084 

-0.00077 

10- 10 

3.58 x 10" 4 

0.961 

0.0084 

-0.00077 


Table 1: Inflationary parameters for Nf = 1, N c = 10, l/g^R — 1 /g\ R = 1, /x = 1 and N = 50. 

Next we discuss the dependence on the renormalized four-fermion coupling, g^ R and the renormal¬ 
ization scale, /x. Since a realistic value for the density fluctuation is obtained for a ~ O(10 -12 ), we set 
a = 10 -12 and numerically calculate the inflationary parameters as g and g^ R vary. We show some 
typical values in Table [2] The /x and g 4 R dependences are also negligible for the density fluctuation, S , 
the spectral index, n s , the tensor-to-scalar ratio, r, and the running of the spectral index, a s . 


5 

l/g4R-l/gt R 

5 

n s 

r 


10~ 8 

1 

3.59 x 10- 5 

0.961 

0.0083 

-0.00076 

10~ 4 

1 

3.59 x 10" 5 

0.961 

0.0083 

-0.00076 

1 

1 

3.58 x 10~ 5 

0.961 

0.0084 

-0.00077 

1 

10“ 4 

3.59 x 10" 5 

0.961 

0.0083 

-0.00076 

1 

10 -8 

3.59 x 10" 5 

0.961 

0.0083 

-0.00076 


Table 2: Inflationary parameters for Nf = 1, N c = 10, a = 10 12 and N = 50. 


It is found that the spectral index, the tensor-to-scalar ratio and the running of the spectral index 
are almost fixed at n s = 0.961, r = 0.0083, a s = —0.00076 for N = 50 in Nf = 1, N c = 10 case. These 
values are consistent with the Planck 2015 data. 
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4.2 Large number of fermion species 

Let us investigate the inflationary parameters at fixed gauge coupling limit, b —> +0 for simplicity. The 
assumption is valid for 11-/V C ~ 2 Nf. However, we vary the number of the fermion species, Nf and N c , 
with a fixed gauge coupling, g. Since the extremely small gauge coupling, a ~ 0(1O~ 12 ), is considered 
the running of the gauge coupling in Eq. (EH) can be neglected. The effective potential develops a non¬ 
trivial extremum for NfN c > 24-7T 2 at the small gauge coupling limit. It induces the fine-tuning problem 
discussed in App. [A] In order to avoid the problem we consider the limited case NfN c < 247r 2 . 



Figure 3: Behavior of the density fluctuation S as 
a function of g for Nf = 1, a = 10~ 12 , l/g 4 R — 

i/4 = i- 
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0.00001 


N f - 1 


N = 60 
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^=*10 


N f = 23 


0.2 0.4 0.6 0.8 1.0 




Figure 4: Behavior of the density fluctuation S as 
a function of /i for N c = 10, a = 10 -12 , l/g^R — 

Wr = 1 - 


In Figs. [3] and [I] we compare the behavior of the density fluctuation between a small and a large 
number of fermion species. The density fluctuation grows up near g = 1 for {Nf,N c } = {1,230} and 
{10,23} cases. Non-negligible g dependence appears for a larger Nf or N c cases. As is shown in Fig. [3] 
the obtained density fluctuation for N c = 230 approaches to the one for N c = 10 as the renormalization 
scale g decreases. We note that lines for N c = 10 and N c = 230 cross twice in g, £ [0,1] and a small but a 
finite difference remains at g —> 0 limit. In Fig. [I] lines for Nf = 1 and Nf = 23 do not cross in g £ [0,1] 
and a larger difference observed at the limit, g —> 0. 



-0.00055 
-0.00060 
^ -0.00065 
-0.00070 



n. 


0.962 0.964 0.966 0.968 0.970 
n s 


Figure 5: Behavior of the tensor-to-scalar ratio r and the running of the spectral index a s as a function 
of the spectral index n s for a = 10~ 12 , Nf = 1, N c = 230 and l/gm — 1 /g% R = 1. 


For a larger number of fermion species the g dependence appears in the spectral index, n s , the tensor- 
to-scalar ratio, r, and the running of the spectral index, a s in models. We draw the behavior of these 
inflationary parameters on the interval g = [0.1,1.0] in Figs. [5] and [6] It should be noticed that we 
normalize the mass scale by the Planck scale and consider the large A limit. Such assumption is suitable 



























-0.00050 
-0.00055 

-0.00060 

>3 

^ -0.00065 
-0.00070 
-0.00075 

0.962 0.964 0.966 0.968 0.970 0.962 0.964 0.966 0.968 0.970 

n s n s 

Figure 6: Behavior of the tensor-to-scalar ratio r and the running of the spectral index a s as a function 
of the spectral index n s for a = 10” 12 , Nf = 10, N c = 23 and l/g±R — 1 /g% R = 1. 

for /i -C 1.0 ~ M p . A similar behavior is observed for {Nf,N c } = {1,230} and {10,23} cases. Thus n„, 
r and a s seem to depend on the product NfN c . At the limit, g —> 0, the result approaches to the one 
for a small number of fermion species, Nf = 1 and N c = 23. The obtained trajectories are shifted as the 
e-folding number varies. All the trajectories are consistent with the Planck 2015 observations. 





Figure 7: Behavior of the density fluctuation 6 as 
a function of N c for Nf = 1, a = 10 -12 , 1/gAR — 
1 / gl R = 1 and g = 1. 



Figure 8: Behavior of the density fluctuation 6 as 
a function of Nf for N c = 10, a = 10 -12 , l/g^R — 

l /gl R = 1 and M = !■ 


The inflationary parameters depend on the renormalization scale g for a larger number of the fermion 
species. It is observed that the parameters, n s , r and a s simply approach to the values for a small number 
of the fermion species at the small g limit. Next we fix the model parameters at g = 1 , a = 10 -12 , 
— 1/s1r = 1 an d evaluate the Nf and N c dependences. The behavior of the density fluctuation, 
S , is plotted as N c and Nf vary in Figs. [T] and [H respectively. The range of the density fluctuation is 
smaller than the one in Fig. [2] The number of the fermion species induces a sub-dominant contribution 
to the density fluctuation compared with the gauge coupling. Thus, we can obtain a suitable density 
fluctuation by tuning the strength of the gauge coupling. 

In Figs. [9] and [10] the tensor-to-scalar ratio, r, and the running of the spectral index, a s , are plotted 
as functions of the spectral index, n s , for the e-folding number, N = 50 and 60. Both the figures show 
a similar behavior to Figs. [5] and [6] A model with a large number of fermion species at g ~ 1 and a 
small number of fermion species at a small g generate nearly equal values of n s , r and a s . The obtained 
trajectories are consistent with the Planck 2015 results. The tensor-to-scalar ratio, r, monotonically 
decreases as N c and Nf increases. The spectral index, n s and the running of the spectral index, a s have 
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n s n s 

Figure 9: Behavior of the tensor-to-scalar ratio r and the running of the spectral index a s as a function 
of the spectral index n s for a = 10” 12 , Nf = 1, 1/gm — 1 / g* iR = 1 and ji = 1. 




n s n s 

Figure 10 : Behavior of the tensor-to-scalar ratio r and the running of the spectral index a s as a function 
of the spectral index n s for a = 10 ~ 12 , N c = 10 , l/gw ~ l/fl , 4 J R = 1 and fJ> = 1 - 


maximal values at a = 2NfN c ~ 380. 

Thus, we demonstrated that realistic gauged NJL-model under discussion for specific values of pa¬ 
rameters is consistent with Planck cosmological data. From another side, the model under discussion is 
realistic high energy physics model of elementary particles. For numerical calculations, we considered 
several choices of fermion species so that such choice is quite standard in particle physics. Thus, we arrive 
to realistic (composite) particle physics model which predicts viable inflation. 

5 Exit from inflation 

The end of the slow roll era is found by observing the parameters, e and tj. In the present model the 
parameter e exceeds unity at er = a enc i in Fig. [lj Under the slow roll approximation, the deceleration 
parameter, q, is written as 


q = 


a 

Wa 


e-1. 


(5.1) 


The deceleration parameter changes the sign at a = a en d and develops a positive value. Thus the universe 
turns to decelerated expansion era. 

In this paper we set the initial condition ctjv to generate the e-folding number between 50 and 60. 
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Thus the duration of the inflation can be estimated by 


tduration — / ( “ 

" &end 


fdaV 


\dtp J dVE/da 


da. 


We numerically find 

UuraUon = (7 ~ 8) X 10 *M~ l ~ 10" 36 S. 
From Eqs. and dsm we obtain 

H 2 = ±(-H + V e ). 

We define the perturbation AH ( t ) by 

H = H dS + A Hit), 

with Has = Ve/ 3 . Inserting Eq. (15.511 into Eq. (15.41) . we get 

A Pf 

— ~-6H dS + 0(AH 2 ). 

The solution of this equation is given by 


A H(t) = A sinh [ 6\j — t 


(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 


where we set an initial condition A H(t = 0) = 0 and A is a constant parameter. Thus the de Sitter 
solution (13.111) is unstable. As is shown in Fig. [1] the potential energy extremely decreases at a = 
o'end- Then the dominant contribution to derive the expansion of the universe comes from the radiation. 
Therefore, the standard cosmology is recovered. It should be noted that R 2 terms induced by the trace 
anomaly also unstabilizes the de Sitter solution and then the radiation dominant era realizes m- 


6 Conclusions 

In summary, we proposed gauged NJL inflation as an alternative for Higgs inflation. Specifically, the 
gauged NJL model is studied as the realistic particle physics model for the composite scalar. Applying 
the auxiliary field method and performing the renormalization group improvement (for review of RG in 
curved spacetime, see m), the model is represented by the gauge-Higgs-Yukawa theory. We assumed 
that the gauge-Higgs-Yukawa theory well describes the model at the scale of inflation at the A —> oo limit. 
Then the produced CMB fluctuations have been calculated in the gauge-Higgs-Yukawa theory with the 
compositeness conditions at the fixed gauge coupling. 

The behavior for inflationary parameters <5, n s , r and a s has been numerically evaluated as a function 
of the couplings, the renormalization scale, the e-folding number and the number of the fermion species. 
It is found that the smallness of the density fluctuation can be realized by tuning the gauge coupling. 
We proved that the model with a small number of fermion species, Nf = 1, N c = 10 satisfies the 
Planck 2015 data. In this case the spectral index, the tensor-to-scalar ratio and the running of the 
spectral index develop almost fixed values, n s = 0.961, r = 0.0083, a s = —0.00076 for N = 50. For a 
large number of fermion species, these inflationary parameters depend on the renormalization scale /i. 
However, the obtained trajectories of the inflationary parameters are consistent with the Planck 2015 
data. Therefore, the gauged NJL model naturally generates inflationary parameters consistent with the 
Planck 2015 observation and may serve as viable inflationary theory. 

At the next step, one has to consider the end of inflation and graceful exit. It is clear that higher- 
derivative RG corrected terms (like R 2 , etc) which are induced as vacuum polarization beyond the 
linear-curvature approximation may naturally provide the graceful exit. It is also interesting to consider 
the reheating process in the gauged NJL model. In fact, this goes beyond this work. However, the 
preliminary considerations of preheating process indicate reheating maybe successful as it is in usual 
scalar theories of inflation (in fact, our model is effectively scalar inflationary theory). These questions 
will be considered elsewhere. Furthermore, NJL model may pretend to the role of dark matter (see, for 
instance, Ref. |22j ). hence presumably unifying inflation with dark matter. 
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A NJL limit 


Let us discuss the slow-roll solution in the NJL model. It is obtained by omitting the gauge interaction. 
Here we consider the NJL limit, a —> 0, where the gauge interaction is eliminated. At the limit the 
running couplings ( 12.011 and ( 12.1011 reduce to 


y%(t) = 


167T 2 a 
2a a c 



A a (t) 

vi(t) 


2 a a c 
167r 2 a 



16tt 2 2 

2a ln(p, 2 /A 2 )’ 


From Eqs. (12.1411 and (IA.Hl we obtain 


m 2 {t) 




(A.l) 

(A.2) 


(A.3) 


where we renormalize 34 (A) suitably to obtain a finite expression |16j . 

Following the similar procedure in Sec. [2] we calculate the RG invariant effective potential at the limit 
A —> 00 and obtain 


1 


1 


V(a) = --- ^ + -a 2 . 


2 V 9iR glu 


R 


12 


Thus, the gravitational effective action with the weak curvature R <C cr 2 becomes 


S (NJL) = j d ± Xy fZT g 


-l R+ lg^d^a-V(a) 


(A.4) 


(A.5) 


By the Weyl transformation m , one can eliminate the interaction term between the curvature and the 
auxiliary field cr. In the Einstein frame the action is simplified as 


S 


(NJL) 


= / d^xsj-^g 


-\r + - Ve(ct) 


(A. 6 ) 


with the effective potential 


V e (<j) 


1 

g4R 


1 

g*4R. 


2 2 
[1 a 


(1 + cr 2 /6) 2 


(A.7) 


The denominator (1 + cr 2 / 6) 2 comes from the NJL limit of f l 2 in Eq. (13.211 . It suppresses the effective 
potential for a large a. The effective potential E3 lias a maximum at cr = \/6. 

In Fig.[lT]the typical behavior of the effective potential (IA.7I) is shown. It is observed that the horizon 
crossing cr = crjy should take place near the maximum of the potential to obtain large enough e-folding 
N = 50 ~ 60. The e-folding is plotted as a function of ctn in Fig. [12] It blows up at cr = \/6. We noted 
that the 34 # dependence gives only a negligible contribution to the e-folding. The result introduces a 
fine- tuning problem why gn ~ \f&- To avoid this fine-tuning problem one should modify the potential. 
Since the potential can be easily modified by including terms which are neglected in our approximation, 
the fine-tuning problem seems to be not so serious. 
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<J 

Figure 11: Behavior of the potential Ve(ct) as a 
function of cr for /i = 1 and 1/ g^R — 1 / g%R = 0.5 x 
10 _15 ,1 x 10 -15 ,2 x 10 -15 . It should be noticed 
that we do not adopt any model parameters whose 
effective potential has a local maximum except for 
this appendix. 



Figure 12: Behavior of the e-folding number N as 
a function of cx/v- 
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